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Abstract 

We derive expansions of the resolvent Rn{x, y; t) = {Qn{x; t)Pn{y; t) — Qn{y; t)Pn{x; t))/{x — 
^ . y) of the Hermite kernel Kn at the edge of the spectrum of the finite n Gaussian Unitary 

^ . Ensemble (GUE„) and the finite n expansion of Qn{x;t) and Pn{x;t). Using these large n 

I expansions, we give another proof of the derivation of an Edgeworth type theorem for the 

'sj" [ largest eigenvalue distribution function of GUE„ . We conclude with a brief discussion on the 

derivation of the probability distribution function of the corresponding largest eigenvalue in 

■ the Gaussian Orthogonal Ensemble (GOE„) and Gaussian Symplectic Ensembles (GSE„). 

o 

. J . 1 Introduction 

■ The author stressed in [1] the importance of having a large n-expansion of the dis- 
tribution of the largest eigenvalue from classical Random Matrix Ensembles. In this 
paper we present another derivation of the probability distribution function of the 
largest eigenvalue from the GUE„. Unlike the previous derivation which follows 
from the Fredholm determinant representation F{XMax < ^) = det(/ — Kn)(t,oo), 
this one follows directly from the resolvent kernel representation F{XMax < = 
exp{— Rn{x, X] t) dx}. (A proof of this representation can be found in [19].) Here 
the Fredholm determinant expansion is replaced by the large n-expansion of P„ and 
Qn- In doing this we discover new integrals relating Painleve functions appearing in 
the study of the largest eigenvalue in Gaussian Ensembles. These large 77,-expansions 
can be used for the analogous problem of finding the probability distribution of the 
largest eigenvalue in the GOE„ and the GSE„ case. 
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Recall that for Gaussian Ensembles, the probability density that the eigenvalues 
are in infinitesimal intervals about the points Xi, ■ ■ ■ , Xn is given by 

FnpiXi,- ■ ■ ,Xn) = C„^exp I j JJ - Xfcl^, (1.1) 

\ 1 / j<k 

with 

— oo < Aj < oo, for z = 1, ■ • • , n, (1.2) 

and Cni3 is the normalization constant. 
Let 

F„,^(t) = P(AL, < t) (1.3) 

be the probability distribution function of the largest eigenvalue in GOE„ for P = 1, 
GUE„ for [3 = 2, and GSE„ for /3 = 4 respectively. 

For the Gaussian n Ensemble, the expected value of the largest eigenvalue is asymptot- 
ically \/2n. Therefore as the size of the matrices grows, so does the largest eigenvalue. 
To have a nontrivial limit, we must center and normalize A^^^,^. In doing this we keep 
the fine tuning constant c introduced in [1] , 

A/3 ._ ALx-(2(r^ + c))^/^ 

max • 2^^/2/1"^/^ ' \ ■ I 

To state our results we need the following definitions. Recall that if 

1 2/ 

' (2'^n!v^)V2 ' 
with Hn{x) the Hermite polynomials of degree n, then the Hermite kernel is 



k=0 

The resolvent of the integral operator on L^(t, oo) with Hermite kernel will be denoted 
by Rn and its kernel denoted by 

Rn{x,y):=iI-Kn)-'Knix,y). (1.5) 

This resolvent also has the representation (see for example [19], page 6) 

Rn[x,y;t) = (1.6) 

x-y 

where 

g„,,(x;t) = iiI-Kn)-\x'ipn) (1.7) 

and 

Pn,i{x; t) = {{I- KnY^ , xVn-l). (1.8) 
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We introduce the following quantities 

qn,ii^) = Qn,i{t;t), Pn,i{t) = PnA^^t) 
Vn,i{i) = iQn,i^^n-l), and Wn,i{t) = {Pn,i,^n-l)- 



(1.9) 

;i.io) 
;i.ii) 



Here (■,■) denotes the inner product on L^(t, oo). In our notation, the subscript 
without the n represents the scaled limit of that quantity when n goes to infinity, and 
we dropped the second subscript i when it is zero. 

If Ai is Airy function, the kernel Kn{x, y) then scaled to the Airy kernel 



KAiiX,Y) 
Our convention is as follow; 



Ai{X) Ai'(y) - Ai(y) Ai'(X) 
X-Y 



Q,(x;s) = ((/-KAi)-\ x^Ai) 
P,{x;s) = ((/-AAi)~\ x^Ai') 



(1.12) 

(1.13) 

(1.14) 
(1.15) 
(1.16) 
(1.17) 



qi{s) = Qi{s;s), pi{s) = Pi{s;s), 

Ui{s) = (Qi,Ai), v,{s) = (P„Ai), 

Vi{s) = {Qi,Ai'), and Wi{t) = {Pi, Ai'). 

Here ( ■ , ■ ) denotes the inner product on L^(s, oo) and ? = 0, 1, 2. 
We use the subscript n for unsealed quantities only. 

Our first result are large n-expansions of Rn{x,y;t), Qn{x;t) = Qn,o{x]t) and 
P„,(x; t) = Pn,o{x; t). 

Theorem 1.1. For 

Y g 

X = y/2{n + c) + y = ^/2{n + c) + -^-^ and t = ^/2{n + c) + ^^, (1.18) 

22n6 22^6 22n6 



as n ^ (yo with X , Y , and s bounded, 
Rn{x,y]t)dx ■■ 



_2 r 

R{X,Y;s)-cQ{X;s)Q{Y;s)n--^ +^ 



Pi{X-s)P{Y- s) + P{X■s)P^{Y■s) 



-Q2{X; s)Q{Y; s) - Qi{X; s)Q^{Y; s) - Q{X; s)Q2{Y; s) + 20c\o{s)Q{X; s)Q{Y; s) 
3 - 20c2 



+ 



{P{X;s)Q{Y; s) + Q{X; s)P{Y; s) 



+ 0(n-i)e„(X,y) 



dX. (1.19) 



The error term, e„(X, Y), is the kernel of an integral operator on L'^{s, oo) which is 
trace class. 



^For the precise definition of this scaUng, see the next section 
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And 

Theorem 1.2. For 



X 



X 



V2(n + c) + -^-^ and t = ^/2{n + c) + 

22?7,6 2277.6 



as n ^ (yo with X and s bounded, 



Qn{x; t) = n<i 



Q{X- s) 



2c -I 



P{X-s)-cQ{X-s)u{s) 



n 3 



+ 



(lOc' - 10c + ^)Qi{X- s) + P2{X; s) + (-SOc' + 10c + ^)Q{X; s)v{s) 



+Pi{X- s)v{s) + P{X- s)vi{s) - Q2{X- s)u{s) - Qi{X- s)ui{s) - Q{X; s)u2{s) 



+ (-lOc^ + ^)P{X; s)u{s) + 20c^Q{X; s)u\s) 



n 3 



20 



0{n~')E^{X;s) 



;i.20) 



and 



Pn{x; t) = nt 



Q{X; s) + 



2c +1 



P{X;s) - cQ{X;s)u{s) 



n 3 



+ 



(10c' + 10c + ^)Qi{X; s) + P2{X- s) + (-30c' - 10c + '^)Q{X- s)v{s) 



-Pi(X; s)v{s) + P{X- s)vr{s) - Q2{X- s)u{s) - Qi(X; s)ui{s) - Q{X; 5)^2(5) 



-10c' + ^)P{X; s)u{s) + 20c'Q(X; s)u\s) 



_ 2 
n 3 



0{n-')Ep{X;s) 



(1.21) 



Theorem II .21 together with the work of Tracy and Widom in [19], all gives another 
proof of the following large n-expansion of (11. 3p when j3 = 2. 



Theorem 1.3. We set 

t = (2(n + c))2 + 2 272 es and 

Ec,2{s) = 2wi - 3u2 + (-20c' + 3)t;o + uiVq - uqVi + uqVq - uIwq. 
Then as n ^ 00 

1 



Fn,2{t) = F2{s)^l + cuo{s)n-^ 
uniformly in s, and 

F2{s) = lim Fn,2it) = exp 



20 



{x — s)q{xY dx 



;i.22) 
;i.23) 

;i.24) 
;i.25) 



is the Tracy- Widom distribution. 

In §2 we derive Theorem 11.11 and Theorem II. 2[ In §3 we give another proof of 
f ll.24p where this time we make use of the representation of the probability distribution 
function of the largest eigenvalue in term of the resolvent of the Hermite kernel instead 
of the Fredholm determinant representation using the Hermite kernel. The advantage 
of this derivation is that this technique also applies to the finite n GOE„ and GSE„. 
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2 Large n- Expansion of Rn{x, y; t), Qn{x] t) and Pn{x] t) 

In this section we will make use of the following expansion of the Hermite kernel 
Kn,2{x,y) derived in [1]. Let Ai{x) be the Airy function and 

, , Ai(x)Ai'(y) - Ai(y)Ai'(x) p,., 

KAi(x,y) = ^-^ — —= / Ai{x + z) Ai{y + z) dz 2.1) 

x-y Jo 

the Airy kernel. For x and y defined by fll.lSI) . we have as n ^ oo 
K„,2(x, y) dx = { AAi(X, Y) - cg Ai{X) Ai{Y)n--^ + 

— [{X + Y) Ai\X) Ai'(F) - (X^ + XY + Y^) Ai{X) Ai(Y)+ 
20 



^ ^ (Ai'(X) Ai(y) + Ai(X) Ai'(F)) 



n-'i + 0{n"^)E{X, Y) | dX, (2.2) 



uniformly in s. The error term, E{X, Y), is the kernel of a trace class integral operator 
on L^(s, oo). 

In order to simplify the notations in this paper, we treat each term appearing 
in (12. 2p as an integral operator as well as the kernel of that operator. For example 
Ai(a;) Ai{y) will be the integral operator with this kernel. 

We recall that the resolvent operator is Rn{x,y]t) = {I — Kn)~^ Kn{x,y). We 
have a large n-expansion of the Hermite kernel Kn{x,y), therefore to derive an ex- 
pansion for the resolvent kernel we only need to derive an expansion of the kernel 
(/ — Kn )~^ {x,y) and multiply the two operators to have our desired result. The 
first part of this section will be devoted to doing that, in the second part we will use 
that result to derive the expansion for The third part will derive an expansion of 
Qn and P„. 

2.1 Large n- Expansion of (/ — )^ {x, y) 

Let 

t{x) = y^2{n + c) + x2-^n-^ (2.3) 

be the scaling function and x the characteristic function of the set {t, oo). 

If L is an integral operator with kernel L{x, y), we will write Lr for the scaled integral 

operator with kernel 

^ rL{r{x),r{y)). (2.4) 



\/2n(i 

With this convention the scaled Hermite kernel (12. 2p has the following representation 
K4T{X),r{y))dT{x) = r' K,,{T{X),r{y)) = KAiiX,Y) - cAi{X) Ai{Y)n~'s + 

— (X + Y) Ai\X) Ai'(y) - {X^ + XY + y2) Ai{X) Ai{Y)+ 
20 
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-20c2 + 3,,., 



(Ai'(X)Ai(y) + Ai(X)Ai'(F)) 



n-^^ +0{n'^)E{X,Y). (2.5) 



Following Tracy and Widom, we denote the kernel of ( / — A'„^2)^^ by pn{x,y), the 
characteristic function of the set (t, oo) by x(^); and the scaled function Xr the 
characteristic function os the set (s, oo) as r(s) = t. Note that 



I - K 



n,2) 



I - Xr Kn,2T Xr) ^ ■ 



(2.6) 



But Xt scales to the characteristic function of the set (s, oo) with t = t{s). To simplify 
notations we will not mention explicitly X{s,oo), but think of the various operators as 
acting on the set (s, oo). With this in mind, we see that the kernel of ( / — Kn,2)~^ 
is ^ 

Pn{x,y) = ("/ - J-^K„,2(r(X),r(F))' 



We combine this with (12. 5p to have 

Pn{x,y) = - KM{X,Y) + cAi{X)Ai{Y)n 

-20c2 + 3 



20 



(X + Y) Ai'(X) Ai'(y) 



■{X^ + XY + Y') Ai{X) Ai(y) + 



(Ai'(x)Ai(y) + Ai(x)Ai'(r)) 



_ 2 

n 3 



+ Oin-')E{X, Y) 



I - KAiiX, Y))-\l + {I - Km{X, Y))-' 



1 

20 



(2.7) 



cAi(X) Ai(y)n" 



{X + Y) Ai'iX) Ai'(F) - {X^ + XY + Y^) Ai{X) Ai{Y)+ 



(Ai'(X)Ai(F)+Ai(X)Ai'(y)) 



n-i + 0(n-^)E(X, Y) 



-1 



(2.8) 



We now think of each term in the large bracket as kernel of an integral operator on 
(s, oo). We know the existence of ( / — Kji^i)~^ . If we factor out this operator in the 
last equation, we find that (12. 8p can be represented by 



/ +(/ - AAi)^'(X,y)|cAi(X)Ai(y)n-i-^ 



(X + y)Ai'(X) Ai'(F) 



-{X^ + XY + Y^) Ai(X) Ai(y) 



-20c2 + 3 



(Ai'(X)Ai(F)+Ai(X) Ai'(y)) 



_ 2 

n 3 



+ 0{n-')E{X,Y)jj ■(/ - KArHx,Y). 
We have the following results (see for example |22j). 



(2.9) 
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If M denotes multiplication by the independent variable, then 

the integral operator M* Ai Ai has kernel X* Ai(X) AiCK), 

the integral operator M* Ai' Ai has kernel X' Ai'(X) Ai(F), and 
the integral operator M* Ai Ai' has kernel X* Ai(X) Y^ Ai' {Y) . 
If we denote by 

p{X,Y) the kernel of the operator (/ — /<Ai)~^ 

on (s , oo) , then using representation fll.l3p we find that the kernel of 

( / - KAi)^^ ■ Ai Ai (2.10) 

(the dot here represent operator multiplication) is 

(p(X,Z) , Z^Ai(Z))(^^^)F^ Ai(y) = Q,{X-s)Y^M{Y), (2.11) 

the kernel of 

( / - i^Ai)"^ ■ M' Ai' Ai (2.12) 

is 

(p(x,z) , z^'Ai'(z))^^^^^y^ Ai(y) = p,(X;.)f^- Ai(y), (2.13) 

and the kernel of 

( / - KmY^ ■ M' Ai' Ai' (2.14) 

is given by 

(p(x,z) , z^Ai'(z))(^^^)y^ Ai'(y) = p,(X;.)r^ Ai'(y). (2.15) 

If we substitute these results in (12.91) . we have 

Pn{x,y) = {^I -|-cQ(X;s)Ai(y)n-U ^ Pi(X; s) Ai'(F) + P(X; s)F Ai'(y) 
-Q2{X- s) Ai(F) - Qi(X; s)F Ai(F) - Q{X- s)Y^ Ai(F) + ^^P(X; s) Ai(y) 



2 



+ + ^ g(X;.)Ai-(y) 



n"§+0(n"i)E(X,y)|^ ■(/ - AAi)~'(X,r). 



(2.16) 

We keep the same notation for the error term which is still a trace class operator 
since the product of the bounded operator (/ — i^Ai)"^ with the trace class operator 
E is trace class. 

Note that the operator n^sL in the braces in fl2.16p is a finite sum of finite rank 
operators, and therefore a trace class operator. We have the following representation 
of our scaled operator 

Kn,2\ = (l - n-^A -(^ - ^Ai] ■ (2.17) 



The trace class limit of the first factor on the right is the identity operator which 
is invertible. Then for large n we can assume that (/ — n~sL) is also invertible. 
This operator therefore admits a convergent (in trace class norm) Neumann series 
expansion for large n of the form 

)~1 oo 
= ^n-H^ = I + n~^L + n^^L^ + 0{n~^)E{X,Y). (2.18) 
fc=0 

We need to find a large ri-expansion of n~iL^. 

1 



-cQ{X;s) Ai(F)n-3 + 



20 



Pi(X; s) Ai'(y) + P{X; s)y Ai'(F) 



-Q2{X; s) Ai{Y) - Qi{X- s)Y Ai{Y) - Q{X- s)Y^ Ai{Y) 
-20c2 + 3 



-20c2 



-PiX;s)Ai{Y) 



+ 



-Q(X;.)Ai'(F) 



n-r^ +0{n~^)E{X,Y) 



(2.19) 



If we use the representation (11.101) . we find that this square is 

{-cQ{X; s) Ai(F)n~^) ■ {-cQ{X; s) Ai{Y)n--^) + 0{n"^)Ei{X, Y) 
= c^QiX; s) ( Ai(Z) , Q(Z, s) ) ^ Ai(y ) n"§ + 0{n~')E, (X, Y) 

= c^Q{X; s) u{s) Ai{Y) n-l + 0{n~^)Ei{X, Y). 
If we substitute ( ]2.20p in the series expansion ( 12.18p . we find that ( 12.160 becomes 



(2.20) 



T'pr.{T{X),T{Y)) = ( /-cQ(X;s)Ai(y)n-3 + - 



Pi{X; s) Ai'{Y)+P{X; s)YAi'{Y) 



-Q2iX; s) Ai(F) - Q,{X; s)Y Ai{Y) - Q{X; s)Y^ Ai(F) + ^r^PiX; s) Ai{Y) 



-20c2 



-g(X; s) Ai'(y) + 20c^Q(X; s) u{s) Ai{Y) 



2 

n 3 



^ ■( / - KAr'{X,Y). 

(2.21) 

+Oin-')E{X,Y)-{I - i^Ai)"'(X,y). 

The notation used for the error term suggests that at each step of the expansion, we 
add to the existing error term all the 0(?7,~^)-terms and rename the error term by 
E{X, Y). This is our first result which we restate as the following Lemma. 

Lemma 2.1. Let pn be the kernel of the operator (/ — Kn,2)'^ on {t, oo), and r the 
transformation defined by (12.31) . Then as n ^ oo with x = t{X) and y = t{Y), 



Pn{x,y) 



I - cQiX;s)Ai{Y)n-^+ - 



Pi{X; s) Ai'(y) + P{X; s)YAi'{Y) 
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-Q2{X; s) Ai{Y) - Qi{X; s)Y Ai{Y) - Q{X; s)Y^ Ai{Y) + 7r^P{X; s) Ai{Y) 



+ QiX; s) Ai'(r) + 20c^Q{X; s) u{s) Ai{Y) 



n 3 



2 



K^O''{X,Y). 



-0{n~^)E{X,Y) 



(2.22) 



uniformly in s. 

The error term E is the kernel of a trace class operator on (s, oo). Here P{X, s) 
Po{X,s), Q{X,s) = Qo{X,s), Qi, Pi, and Q2 are defined in (11.131) and (11.141) 

2.2 Large n-expansion of Rn{x,y) 



In this section we will combine Lemma [?!T] and (11.61) to derive an expansion of Rn{x, y). 
(II. 6p says that 



T'Rn{T{X),T{Y)) 



p„(r(X),r(Z)) , r'i^„,2(r(Z),r(y)) ) (2.23) 

(s,oo) 



First the action of (/ — Kxi) ^ on (12.51) gives 



R{X,Y)-cQ{X-s)A\{Y)n--^+- 



Pi{X- s) Ai'(F) + P{X- s)YAi{Y) 



-Q2{X-s) Ai{Y)-Qi{X-s)Y Ai{Y)-Q{X-s)Y^Ai{Yy 



-20c2 + 3 



{P{X-s) Ai(F) + g(X;5) Ai'(y)) 



n~i +0{n-^)E{X,Y). (2.24) 



Next the action of the first factor in (12.221) can be computed as follow: 
The identity will reproduce (I2.24p . the n~3 term will contribute 

-cQ(X;s)(Ai(Z), i?(Z,y))n-i + ^Q{X- s){A:i{Z) , Q{X- s)) Ai(y)n-i, (2.25) 

2 

and the n~3 term will contribute 

Pi(X;s) {Ai\Z),R{X,Y))+P{X-s) {Z Ai' [Z] , R{Z,Y))-Q2{X- s) (Ai(Z) , i?(Z, y)) 



1 

20 



+ 



-Qi{X-s) {Z Ai{Z) , R{Z,Y)) - Q{X- s) [Z^ Ai{Z) , R{Z,Y)) + 
(p{X- s) (Ai(Z) , i?(Z, Y)) + Q{X- s) (Ai'(Z) , R{Z, Y)) 



-20c2 + 3 



+ 20c2Q(X; s) u{s) (Ai(Z) , R{Z, Y)) 



n 3 . 



(2.26) 



9 



To evaluate the various inner-product appearing in fl2.25p and (12.261) we will make 
use of the following representation R{X, Y) = p{X, Y) — 6{X — Y). Thus 

{Ai{Z),R{Z,Y)) = - Ai{Y) + Q{Y;s), 

{A\'{Z),R{Z,Y)) = - Ai'{Y) + P{Y;s), 

{ZAi{Z), R{Z,Y)) = -Y Ai{Y) + Q,{Y;s), 

{Z''Ai{Z),R{Z,Y)) = -Y^ AiiY) + Q2{Y;s), 

{ZA\'{Z),R{Z,Y)) = -Y Ai'{Y) + P,{Y;s). 

We substitute these values in fl2.25p and fl2.26p . then add all the contributions from 
([221, '!2^ and iKm to obtain Theorem O 

2.3 Large n- Expansion of Qn{x;t) and Pn{x;t) 

In this section we will use Lemma 12.11 together with the expansion of ifni^) and 
Lfn-iix) derivecH in [T] to give a large n-expansion of (5„(x; t) and Pn{x; t) defined by 
(ll.7p and (11.81) respectively. To obtain v?„(x) from Theorem 1.1 of [1], we need to make 
the substitution c — > c + | and the factor in Theorem 1.1 is now ne. For (fn-i{x), we 
need to make the substitution c c + | and the factor is now ne 
We assume without lost of generalities that n = 2k is even, and for the Laguerre 
polynomial of degree k and order a > — 1 we define 



and 



2) (2"n!v/^)i/2 



fc^(_l)fc22fc(A;!)Lfe^(3:')e""'/' 
(22fc(2A:)!V7F)i/2 



Vn-l[X) 



Hn-lix) e 



2) \2^~\n - ly.^yi^- 



We then have for x = t{X), 

M^) = |Ai(X) + ^^^^ 



ki\ 



\fc-lo2fc-l 



(A; - l)!xL|_i(x2)e-^'/2 



Ai{X)n-^i + 



(22fc-2(2fc-2)!v/?)V2 



;i0c2- 10c+^)XAi(X) 



+ X^A^{X) 



n 3 



20 



0(n"^)Ai(X)| 



and 



^n-l{X) 



716 



i|Ai(X) + i^^Ai'(X)n-i 



•^This is the direct consequence of Theorem 1.1 of [I] 



(2.27) 



;ioc2 + ioc+ -)A:Ai(x) 
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_ 2 

n 3 



0{n-^)Ki{X) 



(2.28) 



Next we apply the scaled operator (/ — Kn) ^ to these functions to have our stated 
result. First 

(J - i^Ai)-' ^n(r(X)) = nl (q(X; s) + s)n-^ + 



and 



(10c2-10c+-)Qi(X;s) + P2{X-s) 



_ 2 

n 3 



0{n-^)Q{X-s) 



(10c'-10c+-)Qi(X;s) + P2(X;s) 



n 3 



20 



0{n~^)Q{X-s) 



Preceding in a similar fashion as in the derivation of we make the first factor in 
the right of the fl2.22l) acts on these last two functions and have Theorem II. 2[ Note 
that the inner products here are of the form 

(Z' Ai(Z), Q(Z,s)) = u,{s) and (Z^ Ai'(Z) , Q(Z, s)) = Vi{s). 



To conclude this section, we give the following consequence of Theorem II. 2 [ If we set 
t = r(s) then as n — oo 



Qniris)) = Qnir{s);T{s)) = ( q{s) + 



2c- 1 



p{s) — cq{s)u{s) 



+ (lOc^ - 10c + ^)gi(s) + p2{s) + (-SOc^ + 10c + ^)q{s)v{s) 
+Pi{s)v{s) + p{s)vi{s) - q2{s)u{s) - qi{s)ui{s) - q{s)u2{s) 



and 



-lOc^ + -)p[s)u{s) + 2{]c^q{s)u^{s) 



n 3 



20 



(2.29) 



2c + 1 



p{s) — cq{s)u{s) 



+ 



(10c2 + 10c+-)gi(s)+p2(s) 



-30c2 - 10c+ -)q{s)v{s) 



+Pi{s)v{s) + p{s)vi{s) - q2is)u{s) - qi{s)ui{s) - q{s)u2is) 



+ (-lOc^ + -)p{s)u{s) + 20c'^q{s)u^{s) 



_ 2 

n 3 
~20' 



+ 0{n-')ep{s) 



(2.30) 
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Uniformly in s. We use the notation 

qi{s) = Qi{s;s), pi{s) = Pi{s;s), eg{s) = Eq{s; s) and ep{s) = Ep{s;s) 

and the subscript n is reserved for functions depending on the size n of the matrices 
in consideration. 

3 Large n- Expansion of Fn^2{t) 

In this section we will use the following Fredholm determinant representation of the 
probability distribution function of the largest eigenvalue Fn,2(^) in the GUE„ case: 

Fndi) = H^Ma. <t) = det(J - K„,2). (3.1) 

We also have the following two equations, the proof witch can be found in [19] 

^ logdet(/ - K„,2) = -Rn{t,t;t), (3.2) 

^Rn{t,t-t) = -2qn{t)pn{t). (3.3) 

ot 

Equation (13.31) gives 

— \ogdet{I - Kn,2) = -2 J qn{x)pn{x)dx 

where we used the boundary conditions (g„p„)(cxD) = 0. 

Integration by parts and another use of the boundary conditions gives 

(/ qn{x)pn{x)dx\dy = -2 / {x-t)qn{x) Pn{x) dx] 

(3.4) 

and hence 
Theorem 3.1. 

Fn,2{t) = det{I - Kn,2) = exp (^-2 {x - t)qn{x) Pn{x) dx^ . (3.5) 

Observe that this is the finite n analogue of 01.250 . 
Now we set t = r(s) and x = t{X), then 

roc roc 

2/ ix-t)q4x)pn{x)dx = / (X - s)—qn{T{X))pn{T{X))dX 
Jt Js na 

With the help of fl2.29p and (12.301) . the integrand in this last equation is now 

{X - s)\qn{T{x))pn{T{x)) = 

12 



f{x;s) = {x — s)< q^{x) + 2c 



p{x)q{x) — u{x)q^{x) 



n -i + 



(20c2 + 3)g(x)gi(x) + 



2p2{x)q{x) + {—60c^ + 3)q'^{x)v{x) +2pi{x)q{x)v{x) +2p{x)q(x)vi{x) — 2q2{x)q{x)u{x) 
—2qi{x)q{x)ui{x) — 2(f'{x)u2{x) + (— 60c^ + ?>)p{x)q{x)u{x) + 60c^g^(x)M^(a;) 



+ (20c2 - (x) 



_ 2 



0{n-^)e{x) 



(3.6) 



or 

where 
and 



2 

f{x;s) = (x — s) ( + a(x)n~3 + + 0(?7,"^)e(x) 

20 



a{x) = 2c[p{x)q{x) — u{x)q^{x)~\ 



(3.7) 



b{x) 



(20c^ + 3)g(x)gi(x) + 2p2{x)q{x) + (— 60c^ + 3)q'^(x)f (x) + 2pi{x)q{x)v{x) 



-2p{x)q{x)vi{x) — 2q2{x)q{x)u{x) — 2qi{x)q{x)ui{x) — 2q {x)u2{x) + 60c q {x)u (x) 



-60c^ + 3)p{x)q{x)u{x) + (20c^ - 5)p^{x) 



{3.t 



We use X instead of X here to simphfy notation since X is just a variable of integration. 
We therefore have 



exp { — {x — s)(f{x)dx I exp \ — {x — s)a{x)dxn 3 j exp ( — /(x — s)h{x)dx 



2 

n 3 



exp ~ s)e{x)dx 0{n ^] 



1— I [x — s)a{x)dxn ^ + 2 



nx — s)a{x)dx n"3 + £'a(s) 0(r;,~^) ]• 



2 

n '3 



1- / (x - s)6(2;)d2; + Eb{s)0{n~^) 
20 



1 — i?e(a^) 0{n ^) J exp I — {x — s)q^{x)dx 



1— j {x — s)a{x)dx n 3 -|- j 10 



X — s)a{x)dx 



X — s)h{x)dx 



_ 2 
n 3 



+i..WO(«-)}exp(-/,°:-.),^,.M.) 
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where Ea{s), Eh{s), and i?e(s) are the reminder when expanding the exponential 
functions and Ep{s) is the collection of all the 0{n~^) terms. The second factor of 
this last equality is known as the Tracy- Widom distribution. The first factor will be 
the focus on the reminder of this paper. First we will find a simplification for the n~3 

2 

term, then a simplification of the n~3 term. The error term is a consequence of our 
asymptotic. 



3.1 The n 3 term 

We will use the displayed formulas on page 6 of [22] to simplify this factor. First note 
that if we integrate by parts and use the boundary conditions on a[x), 

XD POO / POO \ 

X — s)a(x)dx = / j a(x)dx\dy. (3-9) 
We have 



00 / POO \ POO / POO 



a{x)dx\dx = — c / I 2q{x){p{x) — u{x)q{x)) dx \dy 

00 / poo \ POO / POO \ 

2q{x)q' (x) dx jdy = —c I j{q^{x)) dx\dy 



{-q\y))dy = -c u'{y)dy = cu{s) . (3.10) 

J s 

Note that this is the term from [T]. 

2 

3.2 The n 3 term 

We will simplify b{x) in two steps. In the first part we will simplify the expression 
containing the constant c, and in the second step simplify the remaining expression. 
The expression proportional to the constant c is 

POO 1 2 POO 

10 l{x — s)a{x)dx —2Q(? {x — s)(^qqi — Sq'^v + Sq'^u^ — 3pqu + p'^){x)dx (3.11) 

_J s J Js 

Equation (13.101) says that the first term is IOc^m^(s). Equation (2.12) of [12] together 
with our our definition of qi give qi{s) = sq{s) — v{s)q{s) + u{s)p{s). If we substitute 
this expression of qi{s) in the second term of fl3.1ip . then it becomes 

■00 

(x — s) {xq^{x) — Aq^{x)v{x) + 'i(f{x)u^{x) — 2p{x)q{x)u{x) + p^{x)^dx 



{xq^{x) — Aq^{x)v{x) + 3g^(x)M^(x) — 2p{x)q{x)u{x) +p^{x)^dxdy. 

In the following steps we integrate this last expression 

x(^{x) — 4:q'^{x)v{x) + Sq^ {x)u^ (x) — 2p{x)q{x)u{x) + p^{x) 
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= —2u{x)p{x)q{x) + 2u^{x)q'^{x) + u^{x)(f'{x) — 2v{x)(f{x) — 2v{x)q^{x) + xq^{x) 

—2v{x)q^{x) + xq^{x) + u{x)p{x)q{x) — u{x)p{x)q{x) + 
= —2 (j){x) —q{x)u{x)) q{x)u{x)+q'^ (x) [u'^{x)—2v{x))+q{x) [xq{x)—2v{x)q{x)+p{x)u{x)) 
+q{x)[xq{x) — 2v{x)q{x) + p(x)u{x)) +p{x)[p{x) — q{x)u{x)) 
= —2q {x)q{x)u{x) + q'^{x)q'^{x) + q{x)p (x) + p{x)q (x) 
= (— u{x) —q^{x)u (x) + {p{x)q{x)) = (— + {p{x)q{x)) 

= {^-u\x))"-v"{x). (3.12) 

The second integral in (13.111) is therefore, 

20c2?;(s) - lOcV(s). (3.13) 

This last expression is due to the fact that the functions ttj, u[, Vi and v[ are zero at 
infinity. The derivation of the various integrals used for (13.121) can be found in |22j . 
We showed that the term containing the constant c simphfies to 

20c\;(s). (3.14) 

Note that this is the same term derived in pp. 
In similar way we show that 



oo roo 




s Jy 



{—Sqiq + 3q'^v — 3p^ + 3upq) (x) dx dy 

noo 
(— 3xg^ + Gq'^v — 3upq + 3upq — 3p^) [x) dx dy 

noo 
{^—3q[xq — 2(fv + up] — 3p[—uq + p]) (x) dx dy 

poo poo poo 

-3 / [qp + pq){x) dx dy = 3 / {v{x)) dxdy = 3v{s). (3.15) 



>s Jy Js Jy 

Suppose that L is the integral of / subject to the boundary condition L(oo) = 0, then 



{x — s)l{x)dx= / L{x)dx. (3.16) 

J s 

Using the following representation (the derivation of which can be found in [19],) 

{—qiq + q'^v—p'^ + upq^{x) = (^Ui — uv + (x) , 

we find that equation (13.151) reduces to the following integral 

) 

(ui — uv + w){x) dx = —v{s), (3-17) 
15 



or that 



Ui{s) — u{s)v{s) + w{s) = V (s) = —p{s)q{s). 

2 

At this stage of the simphfication the n~3 term is 



(3.18) 



(20c2-3)v(s) 



{x—s) [6qqi+2p2q+2piqv+2pqvi—2q2qu—2qiqui—2q^U2—2p^) (x) dx 

POO 

= {20c^-3)v{s) - / {x-s)h{x)dx. 



If we note that 

h{x) = {—Qu[ — 2v'.2 — 2v[v — 2v'vi + 2^2^ + (^uf) + 2u'u2 + 2w') (x) 
= (— 6ui — 2v2 — 2viv + 2u2U + u\-{- 2w) (x), 

2 

then the n~3 term becomes 

(-6^1 - 2v2 - 2viv + 2u2U + ul + 2w) (x) dx. (3.19) 

This is where we stop our simplification of this term. To fully simplify this result 
to match equation (11.241) . we need to deriv^ new integral similar to equation (13.181) 
from the following representation, 

/oo 
(6ni + 2i;2 + 2viv — 2u2U — u\ — 2w) (x) dx = 2wi — 3u2 + uiVq — uqVi (3.20) 

We find that the large n-expansion of the probability distribution function of the 
largest eigenvalue for the GUE„ case is given by the following; if we set 

s 



t{s) = ^/2{n + c) + ^^ 
22n6 



(3.21) 



then as n — > oo, we have 



_ 1 
n 3 



F„,2(r(s)) =F2(s)|l + cu{s)n--^ + 

roo 

(20c^ - 3) t;(s) + / {-6ui-2v2-2viv + 2u2U + ul + 2w){x)dx 



20 

uniformly is s. 



(3.22) 



4 Conclusion 

Our motivation in this paper was to find large ?7,-expansion of g„ and p„. The impor- 
tance of such large n expansions is that they not only give a direct proof of Theorem 
11.31 (GUE„ case), but they are essential ingredients in extending Theorem 11.31 to the 
GOE.„ and GSE„ cases. We will return to these cases in a subsequent paper. 

■^The derivation was simple for the term containing the constant c since we can trace it out. 



16 



Acknowledgements: The author would hke to thank Professor Craig Tracy 
for the discussions that initiated this work and for the invaluable guidance, and the 
Department of Mathematical Sciences at the University of Alabama in Huntsville. 



References 

L. N. Choup. Edgeworth Expansion of the Largest Eigenvalue Distribution Func- 
tion of CUE and LUE IMRN Volume 2006, ID 61049, Pages 1-33. 

P. Deift. Orthogonal Polynomials and Random Matrices: A Riemann-Hilbert 
Approach. American Mathematical Society. Courant Lecture Notes 3, 2000. 

P. Deift, Universality for mathematical and physical systems. Preprint, 



arXiv:math-ph/0603038, 



M. Dieng and C. A. Tracy. Application of random matrix theory to multivariate 
statistics, preprint, Arxiv:math.PR/0603543f 

W. Feller. An Introduction to Probability Theory and Its Applications, Vol.11. 
Second edition, John Wiley, 1971. 

T. M. Garoni, P. J. Forrester and N. E. Frankel. Asymptotic corrections to the 
eigenvalue density of the GUE and LUE. J. Math. Physics, Vol. 46, Issue 10, 
103301 (2005). 

I. Gohberg, S. Goldberg, and M. A. Kaashoek. Classes of Linear Operators, Vol. 
I, volume 49 of Operator Theory: Advances and Applications. Birkhauser, 1990. 

I. Gohberg, S. Goldberg, and M. A. Kaashoek. Classes of Linear Operators, Vol. 

II, volume 63 of Operator Theory: Advances and Applications. Birkhauser, 1993. 

I. C. Gohberg, M. G. Krem. Introduction to the Theory of Linear Nons elf adjoint 
Operators, volume 18 of Translations of Mathematical Monographs. American 
Mathematical Society, 1969. 

H. Hochstadt. The Functions of Mathematical Physics, volume 23 of Pure and 
Applied Mathematics: A series of texts and Monographs. Wiley-Interscience, 
1971. 

I. M. Johnstone, On the distribution of the largest eigenvalue in principal com- 
ponent analysis, Ann. Stats., 29(2):295-327, 2001. 

P. D. Lax. Functional Analysis. Wiley-Interscience, 2002. 

M. L. Mehta. Random Matrices, Revised and Enlarged Second Edition. Academic 
Press, 1991. 

[14] F. W. J. Olver. Asymptotics and Special Functions Academic Press, New York, 
1974. 



17 



M. Plancherel and W. Rotach. Sur les valeurs asymptotiques des polynomes 
d'Hermite. Comm. Math. Helv. 1 (1929)227-254. 

A. Soshnikov. Universality at the Edge of the Spectrum in Wigner Random 
Matrices. J. Stat. Phys., 108(5-6):1033-1056, 2002. 

G. Szego. Orthogonal Polynomials. American Mathematical Society Colloquium 
Publications Volume 23 

C. A. Tracy and H. Widom. Level-spacing distributions and the Airy kernel. 
Commun. Math. Physics, 159:151-174, 1994. 

C. A. Tracy and H. Widom. Fredholm determinants, differential equations and 
matrix models. Commun. Math. Physics, 163:33-72, 1994. 

C. A. Tracy and H. Widom. On orthogonal and symplectic matrix ensembles. 
Commun. Math. Physics, 177:727-754, 1996. 

C. A. Tracy and H. Widom. Correlation functions, cluster functions, and spacing 
distributions for random matrices. J. Stat. Phys., 92(5-6):809-835, 1998. 

C. A. Tracy and H. Widom. Airy kernel and Painleve II. In Isomonodromic de- 
formations and applications in physics, volume 31 of CRM Proceedings & Lecture 
Notes, pages 85-98. Amer. Math. Soc, Providence, RI, 2002. 

C. A. Tracy and H. Widom. Distribution functions for largest eigenvalues and 
their applications. In Proceedings of the International Congress of Mathemati- 
cians, Beijing 2002, Vol. I, ed. LI Tatsien, Higher Education Press, Beijing, 
pgs. 587-596, 2002. 

C. A. Tracy and H. Widom. Matrix kernels for the Gaussian orthogonal and 
symplectic ensembles. Ann. Inst. Fourier, Crenoble, 55, 2197-2207, 2005. 

E. T. Whittaker and G. N. Watson. A Course of Modern Analysis. Fourth 
Edition Cambridge University Press, 2004. 



18 



